Ultra-fast spin avalanches in crystals of molecular magnets in terms of magnetic 
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Recent experiments (Decelle et al., Phys. Rev. Lett. 102, 027203 (2009), Ref. fl|) discovered 
an ultra-fast regime of spin avalanches in crystals of magnetic magnets, which was three orders of 
magnitude faster than the traditionally studied magnetic deflagration. The new regime has been 
hypothetically identified as magnetic detonation. Here we demonstrate the possibility of magnetic 
detonation in the crystals, as a front consisting of a leading shock and a zone of Zeeman energy 
release. We study the dependence of the magnetic detonation parameters on the applied magnetic 
field. We find that the magnetic detonation speed only slightly exceeds the sound speed in agreement 
with the experimental observations. 
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Molecular magnetism is a rapidly developinff interdis- 
ciplinary research area within material science 0,3- One 
of the widely investigated materials in the subject is 
Mni2-acetate, with a high spin number (S = 10) and 
strong magnetic anisotropy 2-6]. At sufficiently low tem- 
perature all the spins of the molecules become oriented 
along an applied external magnetic field, thus occupying 
the ground state (e.g. Sz = 10); in this state the magne- 
tization reaches its saturation value. When the magnetic 
field direction is switched to the opposite one, the for- 
mer ground state becomes metastable with an increased 
potential energy (the Zeeman energy) and a barrier sep- 
arating it from the new ground state. Active research 
on the subject demonstrated that spin-relaxation from 
the metastable to the ground state often happens in the 
form of a narrow front spreading in a sample with ve- 
locity of a few meters per second 0-13]. Still, all these 
works focused on magnetic deflagration, i.e. a front of 
energy release propagating due to thermal conduction at 
velocities much smaller than the sound speed. 

However, in contrast to other studies, recent experi- 
ments by Decelle et al., Ref. [H, discovered a new fast 
regime of the magnetic avalanches in Mni2-acetate with 
a front velocity estimated to be (2000-3000) m/s, which 
exceeds the typical magnetic deflagration speed by three 
orders of magnitude. Though a limited number of sensors 
led to rather large uncertainty in measuring the front ve- 
locity, the experiments still indicated clearly that it was 
comparable to the sound speed in the crystals. Further- 
more, Decelle et al. Q hypothetically interpreted the 
new regime as magnetic detonation. Although this hy- 
pothesis looked reasonable, it still required much theo- 
retical work to be justified. In particular, detonations in 
combustion problems demonstrate a propagation veloc- 
ity larger than the sound speed by order of magnitude 
and destructively high pressure 14, l^l- Besides, even in 
combustion science, the phenomenon of deflagration-to- 
detonation transition has remained one of the least un- 
derstood processes for more than seventy years, despite 
its extreme importance 3, It is only recently that 



a quantitative theoretical understanding of this process 
has been achieved ^17iri2Q] . 

In this Letter we demonstrate the possibility of mag- 
netic detonation, in the form of a front with a leading 
shock and a zone of Zeeman energy release. We study 
the dependence of the magnetic detonation parameters 
on the applied magnetic field. We find that the magnetic 
detonation speed is only slightly greater than the sound 
speed, in agreement with experimental observations. 

In line with the experiments [l|, we consider magnetic 
avalanches in Mni2-acetate with the Hamiltonian 



(1) 



suggested in [ll|. Here Sz is the spin projection, j3 « 
0.65i^ is the magnetic anisotropy constant, g m 1.94 
is the gyromagnetic factor, is the Bohr magneton, 
and Hz is the external magnetic field. The Hamiltonian 
(1) determines the Zeeman energy release and the en- 
ergy barrier of the spin transition (in temperature units), 
which depend on the magnetic field as 
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respectively, with Sz — 10. The energy barrier decreases 
with the field while the Zeeman energy increases linearly. 
Next, we consider a stationary magnetic detonation in a 
crystal of molecular magnets. Similar to the theory of 
shocks in solids [2l| , we adopt the reference frame of the 
detonation front and find the conservation laws of mass, 
momentum and energy according to 
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where the label designates the initial state, D is the 
detonation speed, u is velocity produced by the deto- 
nation, e is thermal energy per molecule, and a is the 
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fraction of molecules in the metastable state. Here we 
neglect the thermal conduction, since this is a compar- 
atively slow process. In the theory of shock waves one 
introduces the volume per unit mass V = 1/ p instead of 
density. The conservation laws Eqs. (jH)-® have to be 
complemented by an equation of state. Following Ref. 
2l|, we represent the pressure and energy of condensed 
matter at low temperature as a combination of elastic 
and thermal components according to 
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where cq is the sound speed (we take cq ~ 
accordance to [l|), the power exponent n 
gested in [2l|, F f« 2 is the Gruneisen coefScient, 9_d is 
the Debye temperature with Qd — 38if for Mni2, ks is 
the Boltzmann constant, A = 127r*/5 corresponds to the 
simple crystal model, a = 3 is the problem dimension. 
Thus, in Eqs. (H])-® we have a complete system for 
describing magnetic detonation in molecular magnets. 

The properties of shocks and detonations are repre- 
sente d by the Hugoniot/detonation curve P — P{V), see 
Ref. [3, [H, HH . We also introduce the scaled density 
ratio r = p/po = Vq/V which characterizes the matter 
compression. Using Eqs. (H])-®, we derive the following 
implicit form for the detonation relation 
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In the case of zero energy release (a — 1), Eq. ([9|) reduces 
to the Hugoniot equation for a shock wave (which we de- 
note by the subscript s). In the detonation, the leading 
shock compresses the sample, increases temperature and 
hence facilitates the spin reversal with the Zeeman energy 
release. The released Zeeman energy provides expansion 
of the medium, which acts like a piston and supports the 
leading shock. In the case of the completed spin rever- 
sal (a = 0), Eq. ^ describes the final state behind the 
detonation front (which we denote by the subscript d). 
The insert of Fig. 1 shows the Hugoniot and detonation 
curves found using Eq. ([9|) for _ff = 4 T. We assume that 
there is no external atmospheric pressure and the initial 
temperature is negligible, which corresponds to the ini- 
tial point {V = Vq\ r = 1\ P = {)). Because of the energy 
release, the detonation curve is always above the Hugo- 
niot one. In the case of Mni2 we find that the elastic 
contribution to the pressure and energy dominates over 
the thermal one, which leads to a rather weak detonation 
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FIG. 1: The insert: Traditional presentation of the Hugoniot 
and detonation curves and the tangent line to the detonation 
curve in Mni2-acetate for the external magnetic field = 
4T. The main plot: The Hugoniot and detonation curves with 
the tangent line extracted; label "t" stands for tangent. 



with the shock and detonation curves almost coinciding 
as shown at the inset of Fig. 1. A self-supporting detona- 
tion corresponds to the Chapman- Jouguet (CJ) regime, 
for which velocity of the products in the reference frame 
of the front is equal to the local sound speed [13] . The 
CJ point at the detonation curve is determined by the 
tangent line connecting the initial state and the detona- 
tion curve. Since the detonation and Hugoniot curves 
are extremely close at the insert of Fig. 1, the intersec- 
tion of the tangent line cannot be seen in the traditional 
representation of the curves. In order to make the fig- 
ure illustrative, we subtract this tangent line from the 
Hugoniot and detonation curves in Fig. 1. In the new 
representation, the tangent line corresponds to the zero 
line, while the Hugoniot and detonation curves may be 
distinguished quite well. The CJ point in Fig. 1 corre- 
sponds to the final state behind the detonation front. The 
shock point indicates the strength of the leading shock 
as determined by the Zeeman energy release for the CJ 
regime; the density and the pressure acquire maximum 
values at the shock front. The Zeeman energy release 
behind the shock produces expansion of matter with an 
ensuing pressure reduction. 

We notice from Fig. 1 that the Mni2 crystal is com- 
pressed by few percents in the detonation wave, which 
makes an analytical theory for the detonation front pa- 
rameters possible using expansion r = 1 + 5 with 5^1. 
Then, to leading order in 5, Eqs. ([7]), © reduce to 
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We find the final compression behind the detonation front 
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The compression behind the leading shock is larger by a 
factor of 2, Ss ~ 2(5^, as may be seen from the parabolic 
shape of the Hugoniot and detonation curves in Fig. 1. 
The detonation speed may be found from Eq. Q as 
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which means that the magnetic detonation speed slightly 
exceeds the sound speed in agreement with the experi- 
mental observations [1[. Substituting Ss, Sd, and a = 1; 
into Eqs. (ITUl) . ([TT|). we find the analytical formulas for 
pressure and temperature at the shock and behind the 
detonation front, respectively. Taking into account Eq. 
([3]), these formulas specify the dependence of the detona- 
tion parameters on the external magnetic field. In Fig. 2 
we compare the analytical theory to the numerical solu- 
tion to Eq. ([9]) . The curves of Fig. 2a show the density at 
the shock wave and behind the detonation front, when all 
the spins have been aligned along the magnetic field. We 
see that the density at the shock increases by less than 3 
percents if the magnetic field is below 10 T. Because of 
this small compression, the analytical theory ([T2| is in a 
very good agreement with the numerical solution. The 
maximum value of shock pressure is about 1.2 atm for 10 
T. Thus, due to the small compression and the moder- 
ate pressure increase, the magnetic detonation does not 
destroy the magnetic properties of the crystals. 

At the same time, the crystal temperature increases 
considerably because of the shock, and this stimulates 
a fast spin reversal and a further temperature increase. 
Fig. 2b illustrates the temperature increase at the shock 
wave and behind the detonation front. Again, we observe 
very good agreement between the analytical theory and 
the numerical solution. The temperature at the leading 
shock is comparable to that expected for the magnetic 
deflagration |7HlO|. which also makes the reaction time 
comparable in both processes. In classical combustion, 
the temperature at the leading shock in the detonation 
wave is still quite small in comparison with the activation 
energy of the chemical reactions, so that the active reac- 
tion zone lags considerably behind the shock 15|. The 



situation may be quite different in magnetic detonation. 
When the magnetic field is stronger than 2 — 3 T, the 
shock temperature is relatively high {Ea/Tg < 5) so that 
active spin reversal starts right at the shock wave. Figure 
2b presents also the energy barrier as a function of the ex- 
ternal magnetic field. The energy barrier decreases with 
the growth of the magnetic field, Eq. ^ , as shown in Fig. 
2b. When the magnetic field exceeds 10 T, the energy 
barrier vanishes, the metastable state turns unstable, and 
the molecules may settle down freely to the ground state. 
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FIG. 2: Density ratio (a) and temperature (b) at the leading 
shock and behind the detonation front versus the external 
magnetic field. Solid lines show exact numerical solution; the 
dashed lines stand for the analytical theory. 



Hence, one may interpret magnetic avalanches as deto- 
nation or deflagration only for the fields below 10 T. 

Finally, we describe the internal structure of the mag- 
netic detonation front. In the reference frame of the mov- 
ing front, the molecule fraction with the spin opposite to 
the field direction is determined by 11 1 
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dx 
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where th is a constant of time dimension characterizing 
the spin reversal. We integrate Eq. ([T^ numerically to- 
gether with Eqs. ^ and ([7]) along the tangent line in 
Fig. 1, from the shock to the CJ point; the obtained pro- 
files are depicted in Fig. 3 for = 3 T. The background 
shading represents the energy release due to the spin re- 
versal; the temperature and the pressure are scaled to 
their maximal values. The coordinate is scaled by the 
characteristic length Lq — cqtr w 2 • lO^^m, where we 
take Tfi « 10~^s as obtained in several experiments 0- 
[lo| . Using this value we can estimate the characteristic 
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FIG. 3: Stationary profiles of the scaled temperature, pres- 
sure, and fraction of molecules in the metastable state for 
Hz = 3 T. The background shading shows the energy release. 



width of the stationary detonation wave to a few mil- 
limeters. The applied magnetic field influences strongly 
the reaction rate and thus the front width. For mag- 
netic fields higher than 5 T, the detonation width is < 1 
mm, while for a weaker field the width may increase con- 
siderably. For this reason, the detonation mechanism in 
molecular magnets may only be observed in experiments 
utilizing high enough magnetic fields, since the typical 
sample size is of order of several millimeters. The typical 
scales in the experiments of Ref. [H were also about a 
few millimeters. Thus the experimentally observed fast 
avalanche regime was, presumably, a non-stationary det- 
onation in the process of developing. 

To summarize, in this Letter we have developed a the- 
ory of magnetic detonation in molecular magnets, which 
explains a new regime of ultra-fast spin avalanches dis- 
covered recently in the experiments of Ref. [f*] . The det- 
onation regime is two to three orders of magnitude faster 
than the magnetic deflagration observed before 
We have shown that the leading shock triggers the spin 
reversal in these magnetic systems, and that the magnetic 
detonation propagates with velocities slightly larger than 
the sound speed. In contrast to traditional detonations in 
combustion, which are characterized by strongly super- 
sonic velocities and ultra-high pressure, magnetic deto- 
nations involve rather moderate pressure increase, which 
is about I atm even for considerable magnetic fields. For 
this reason, magnetic detonation does not destroy mag- 
netic properties of the crystals, a very important conclu- 
sion in view of possible applications of molecular magnets 
to, e.g., quantum computing. 
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